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ABSTRACT 

Gas giant planets orbiting within 0.1 AU of their host stars, unlikely to have formed in situ, are 
evidence for planetary migration. It is debated whether the typical hot Jupiter smoothly migrated 
inward from its formation location through the proto-planetary disk or was perturbed by another body 
onto a highly eccentric orbit, which tidal dissipation subsequently shrank and circularized during 
close stellar passages. Socrates and collaborators predicted that the latter class of model should 
produce a population of super-eccentric proto-hot Jupiters readily observable by Kepler. We find a 
paucity of such planets in the Kepler sample, disagreeing with the theoretical prediction with 98.7% 
confidence. Observational effects are unlikely to explain this discrepancy. We find that the fraction of 
hot Jupiters with orbital period P > 3 days produced by the stellar binary Kozai mechanism does not 
exceed 0.15^q'^. Our results may indicate that disk migration is the dominant channel for producing 
hot Jupiters with P > 3 days. Alternatively, the typical hot Jupiter may have been perturbed to 
a high eccentricity by interactions with a planetary rather than stellar companion and began tidal 
circularization much interior to the ice line after multiple scatterings. A final alternative is that tidal 
circularization occurs much more rapidly early in the tidal circularization process at high eccentricities 
than later in the process at low eccentricities, contrary to current tidal theories. 
Subject headings: planetary systems 



1. INTRODUCTION 
Roughly 1% of Sun- like stars host hot Jupiters, giant 



planets with small 


semi-major axes (Mayor et al. 


2011 


Howard et al.|2012 


Wright et al.|2012|). Unlikely to have 



lence of planetary migration, which may take place via 
interactions with the proto-planetary disk (e.g. Goldre 
Tremaine|19801|Ward|1997ljAlibert et al "|2005[|lda 



ich 



& Lin||2008| IBromley & Kenyon||2011 ), or other bodies 



in the system. One or more companions can create a hot 
Jupiter by perturbing a cold Jupiter onto an eccentric or- 
bit, which tidal forces shrink and circularize during close 
passages to the star. Proposed mechanisms for this "high 
eccentricity migration" (HEM) include Kozai oscillat ions 
induced by a distant stellar binary companion (e.g. |Wu 



fc Murray|[2003l |Fabrycky fc Tremam^[2007| ; LN aoz ■ et al 

dNaoz et al 



2012 



2011 



et al 



or by another planet in the system 
Lithwick & Naoz 20 ill), planet-planeF scattering 
Rasio & Ford||1996 IFord & Rasiol 20061 iChatterjeel 
2008 1 fr'ord fc Rasfo||2008[ | J uric fc Trem aine 2008; 



Matsumura et al.||2010[ |Nagasawa fc Ida||2011[ |Beauge| 
fc Nes vorny 2012; Bole y et al.||2012[ ), and secular chaos 
( |Wu fc Lithwick||201ip . 

One way to distinguish whether disk migration or HEM 
is dominant in setting the architecture of systems of gi- 
ant planets is to search for additional populations of gi- 
ant planets that may also result from HEM, including 
1) failed hot Jupiters, which are stuck at high eccentric- 
ities but with periapses too large to undergo significant 

1 Harvard-Smithsonian Center for Astrophysics, 60 Garden St, 
MS-10, Cambridge, MA 02138 

2 rdawsonOcfa.harvard.edu 

3 Division of Geological and Planetary Sciences, California In- 
stitute of Technology, 1200 East California Boulevard, MC 170-25, 
Pasadena, CA 91125, USA 

4 NASA Exoplanet Science Institute (NExScI), CIT Mail Code 
100-22, 770 South Wilson Avenue, Pasadena, CA 91125 



tidal circularization over the star's lifetime, 2) Jupiters 
on short-period, moderately-eccentric orbits, nearing the 
end of their HEM journey, and 3) proto-hot Jupiters on 
super-eccentric orbits in the process of HEM. Recently, 



Socrates et al. 



^ 2012[ ) (S12 hereafter) suggested that, if 
HEM is the dominant channel for producing hot Jupiters, 
we should readily detect a number of super-eccentric 
Jupiters in the act of migrating inward. Moreover, they 
showed that the number of super-eccentric Jupiters can 
be estimated from the number of moderately-eccentric 
Jupiters that have similar angular momentum, based on 
their relative circularization rates. Based on the num- 
ber of moderately-eccentric, short-period Jupiters found 
by other planet hunting program s (tabulated in the Exo- 
planet Orbit Database, EOD, by |Wright et al.|2011[ ), S12 
predicted that the Kepler Mission should discover 5-7 
proto-hot Jupiters with eccentricities e > 0.9 and noted 
tha t these planets shoul d in fact already be present in 
the Borucki et al. ( 2011[ ) candidate collection. 

The S12 prediction requires a steady production rate 
of hot Jupiters throughout the Galaxy, as well as sev- 
eral conventional assumptions, including conservation of 
the hot Jupiter's angular momentum, tidal circulariza- 
tion under the constant time lag approximation, and the 
beginning of HEM at or beyond the ice line. This pre- 
diction is a useful, quantitative test for discerning the 
origin of hot Jupiters. Confirmation of their prediction 
would reveal that hot Jupiters are produced by interac- 
tions with companions, not a disk, while a paucity of 
proto-hot Jupiters in the Kepler sample would inform 
us that HEM is not the dominant channel, or that some 
aspect of our current understanding of HEM is incorrect. 

Motivated by the S12 prediction, we have been using 
what we term the "photoeccentric effect" to measure in- 
dividual eccentriciti es of Jupiter-sized planet s from their 



transit lig ht curves (Dawson fc Johnson|2012 DJ12 here- 
Dawson et al.| ( |2012| ) (D12 hereafter) identi- 



after) 



2 



Dawson et al. 



fied KOI-1474.01 as a transiting planet candidate with 
a long orbital period (69.7 days), a large eccentricity 
(e = 0.81 ± 0.10), and transit timing variations caused 
by a massive outer companion. However, uncertainty in 
the candidate's eccentricity made it ambiguous whether 
KOI-1474.01 is one of the proto-hot Jupiters predicted 
by S12 or, alternatively, a failed-hot Jupiter beyond the 
reach of tidal circularization over its host star's lifetime. 

Here we examine the entire sample of Kepler Jupiters 
to assess whether the planets expected from HEM are 
present. We find with 98.7% confidence that the pu- 
tative highly-eccentric progenitors of hot Jupiters are 
partly or entirely missing from the Kepler sample. In 
Section [2j we summarize the S12 prediction and assump- 
tions. In Section [3j we update the S12 prediction, ac- 
counting for counting uncertainties and incompleteness, 
and translate it into a prediction for transit light curve 
observables. In Section [4j we compare the prediction of 
Section [3] to the light curve properties of candidates in 
the Kepler sample and conclude that there is a paucity 
of proto-hot Jupiters. We consider observational causes, 
finding that they are unlikely to explain the discrepancy 
between theory and observations. In Section [5] we place 
an upper-limit on the fraction of hot Jupiters created by 
stellar binaries, consider the contribution of disk migra- 
tion to the hot Jupiter population, and present Monte 
Carlo predictions for other dynamical scenarios, finding 
that the paucity of proto-hot Jupiters can be compati- 
ble with HEM. We conclude (Section |6| by outlining the 
theoretical and observational pathways necessary to dis- 
tinguish the dominant channel for hot Jupiter creation. 

2. PREDICTIONS AND ASSUMPTIONS BY SOCRATES 
AND COLLABORATORS 

S12 predicted that the Kepler mission should discover 
a number of super-eccentric, hot Jupiter progenitors in 
the process of high eccentricity migration (HEM). Pre- 
viously (DJ12), we showed that super-eccentric plan- 
ets should be easily identifiable from their transit light 
curves and thus precise radial-velocity (RV) follow-up is 
not necessary. This is fortunate as most Kepler stars are 
too faint to be amenable to precise RV observations. To 
predict the number of super-eccentric Jupiters, S12 con- 
sidered a population of proto-hot Jupiters undergoing 
tidal circularization along a "track" of constant angular 
momentum. Using an eccentricity-dependent tidal cir- 
cularization rate (described below), they computed the 
number ratio of super-eccentric to moderately-eccentric 
proto-hot Jupiters along the track. In Section |2.1| we 
follow S12 to derive a formula for the expected number 
of super-eccentric proto-hot Jupiters. In Section 2.2 we 
summarize the assumptions on which the S12 prediction 
depends and how these assumptions affect the expected 
number of proto-hot Jupiters. 

2.1. Number of super- eccentric Jupiters along an 
angular momentum track 

To predict the number of super-eccentric proto-hot 
Jupiters, S12 assumed a steady production of hot 
Jupiters and assessed the relative amount of time spent in 
the early, high-eccentricity phase by a proto-hot Jupiter 
undergoing tidal circularization, compared to the time 
spent at moderate eccentricities later in the process. 
They assumed that the planet's specific orbital angular 
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Table of tidal quantities 
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momentum h is conserved in HEM, following a "track" 
defined by a constant ag na i, the semi-major axis the 
planet reaches once its orbit has fully circularized. For a 
planet with mass M p orbiting a star of M*: 



flfinai = h 2 /[G(M+ + M p )] = o(l - e 2 



(1) 



where G is the universal gravitational constant, a the 
instantaneous semi-major axis, and e the instantaneous 
eccentricity. The angular momentum can also be defined 
in terms of the final orbital period Pflnal' 



P 



final 



(2ira: 



final)/* 1 



P(l 



= 2n3/2 



(2) 



where P is the instantaneous orbital period. 

The number of super-eccentric Jupiters (N sup ) along 
a track of constant angular momentum is related to the 
number of moderately-eccentric Jupiters (-/V m od) by: 



Nsup = N mod r(e ri 







(3) 



where r(e max ) is the ratio of time spent at super- 
eccentricities (0.9 < e < e max ) to moderate eccentricities 
(0.2 < e < 0.6). We place bars over -/V sup and N mod to in- 
dicate that these are mean numbers. The observationally 
counted numbers are sampled from Poisson distributions 
defined by these means. The ratio is 
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set by the timespan P max of the survey. For example, for 
-Pmax = 2 years and P fina i = 5 days, e max = 0.982. The 
eccentricity damping rate e due to tides raised on the 
planet under the constant tidal tim e la g approximation 



Socrates & Katz 



(Egglcton et al.|1998j Hansen|2010j S12,_ 
2012| Socrates et al.|2012 ) , assuming the planet's spin is 
pseudo-synchronous, is: 



(1 - e 2 ) 3 / 2 / e 

CphysOfinai 



(5) 



where f e , a function of e of order unity, and C p h ys , a 
constant based on M*, M p , and the planet's radius R p , 
are defined in Table [TJ 

Most Jupiters in the Kepler sample lack measured ec- 
centricities, and therefore N mo d of the Kepler sample 
is unknown. However, S12 made use of a "calibration" 
sample (which we denote with the subscript 0) of Jupiters 
detected via the transit-method by non-Kepler survey s 
(Figure [TJ compiled from the EOD, |Wright et al.|[20T2 ). 
Because they were detected by transits, the calibration 
sample suffers from the same geometric transit biases as 
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the Kepler sample. As noted by S12, the transit prob- 
ability along a given angular momentum track is con- 
stant. To estimate iV mo( j along a track in the Kepler 
sample, we use the ratio of N mo ^o to the number in 
a class of calibration object that can also be counted 
in the Kepler sample. Ideally, this class would along a 
i-finai track. However, because the eccentricities of the 
Kepler planets are unknown, instead the class we use is 
planets with orbital period P = Pfi na i, of which there 
are N p=pfi n ai,o m the calibration sample. If we assume 
this ratio iV moc i j o/-/Vp = pfj na i i o is the same for calibration 
sample as for the Kepler sample, then we can compute 
the expected N mo & for the Kepler sample: 



Amod — 



N 



mod,0 



N P= 



-N P= 



P=Pfinal 



(6) 



P=Pfinal,0 



Based on Equation (13| 5 S12 predicted 5-7 super 



eccentric proto-hot Jupiters in the Borucki et al 



sample. We will recompute this prediction in Section 
incorporating the latest survey samples, Poisson uncer- 
tainties, and incompleteness. 

2.2. Summary of assumptions forming the basis for the 
S12 prediction 

Here we summarize both stated and unstated assump- 
tions of S12 and infer how violations would affect the 
expected number of super-eccentric proto-hot Jupiters 
in the Kepler sample. Certain assumptions, if vio- 
lated, may result in fewer than expected super-eccentric 
Jupiters. We discuss these assumptions in detail in Sec- 
tion 5: hot Jupiters of orbital periods up to 10 d ays 
have migrated via tidal circularization (Section |5.1[ ); a 
proto-hot- Jupiter typically begins its HEM journey at 
or beyond the ice line, after which it experiences no 
perturbations that per man ently change its angular mo- 
mentum (Sections 5.2 and 5.3 1; a steady "current" of 
proto-hot Jupiters is being prod uced around the sample 
of observable stars (Section 5.4); and the ratio used in 
Equation (|6| is t he s ame for the Kepler and calibration 
samples (Sectio n |5.4[ ). 

In Appendix |A[ we describe additional assumptions. 
S12 assumed that the planet's radius does not change, an 
assumption which, if violated, would not result in fewer 
super-eccentric Jupiters. S12 made several assumptions 
that we do not expect to be violated: angular momentum 
is not exchanged between the star and planet or between 
the planet's spin and orbit, moderately-eccentric Jupiters 
in the calibration sample truly have e > 0.2, and the 
Kepler false positive rate is low. 

Finally, S12 assumed that tidal evolution (for tides 
raised on the planet; see Appendix [A] regarding tides 
raised on the star) occur s according to the constant tidal 
time lag approxim ation jHut||198lT|Eggleton et ah]|1998 
Socrates & Katz||2012[ 



[Socrates et al 
sumption controls the ratio r(e 



20l2[ ). This as- 
oF high-eccentricity 
proto-hot Jupiters to moderate-eccentricity hot Jupiters 
along a given angular momentum track. The constant 



5 S12 used Np— pflnal f° r the Kepler sample but, for the calibra- 
tion sample, the total number of planets observed along the P nna i 
track. However, we wish to treat both samples the same and thus 
use the same type of quantity for both: 7Vp = pfl na i for the Kepler 
sample and 7Vp = p nn ai o f° r the calibration sample. 



tidal time lag approximation is conventional but may be 
violated: if the dissipation rate were larger for highly- 
eccentric Jupiters along a given angular momentum track 
than for moderately-eccentric Jupiters, we would expect 
fewer super-eccentric proto-hot Jupiters than predicted 
or vice versa. 

Dynamic tides, in which dissipation occurs through 
surface gravity waves (e.g. Zahn||1975 1, may be impor- 
tant in a pr o to-hot Jupiter's tidal evolution. |Beauge| 
& Nesvorny (2012) argue that dynamic tides act at 
high eccentricities and equilibrium tides at low eccen- 
tricities; they added an empirical correction factor to 
the constant time-lag model so that, at large eccen- 
tricities, it matches the numeric al results of the dy- 
namic al tide model computed by |Ivanov fc Papaloizou| 
(20111. The empirical correction factor is prop ortional 

to 10 20 ° 9e , where q is the periapse distance (Beauge 



fc Nesv orny 2012 
track, q 



Along a given angular momentum 
-e), so the tidal dissipation timescale 



afinal7P 

[proportional to io 200afi n a i e /( 1 + e )] is longer for larger ec- 
centricities. If this correction factor applies, the contri- 
bution of dynamical tides would increase the expected 
number of super-eccentric proto-hot Jupiters. 

The effect of tides on orbital evolution remains uncer- 
tain and is a topic of ongoing research. It remains unclear 
whether the true typical tidal evolution would result in 
more or in fewer super-eccentric Jupiters. Our results 
should be revisited as this subfield continues to advance. 

3. UPDATED PREDICTION FOR NUMBER OF 
SUPER-ECCENTRIC PROTO-HOT JUPITERS AND 
TRANSIT LIGHT CURVE OBSERVABLES 



In Section [3.1| we derive the expected number of identi- 
fiable Kepler super-eccentric proto-hot Jupiters, follow- 
ing S12 but using updated survey samples. We refine 
the S12 prediction by quantifying its uncertainty and in- 
corporating incompleteness due to the limited timespan 
of the data. In Section 13.21 we describe how to con- 



firm or rule out the existence of super-eccentric proto-hot 
Jupiters using Kepler photometry alone by recasting the 
prediction in terms of light curve observables. 

3.1. Expected number of proto-hot Jupiters with e > 0.9 
in the Kepler sample 



In Section |2.1[ we followed S12 to derive an equa- 
tion for the expected number of super-eccentric proto-hot 
Jupiters for a given Pfi na i (Equation [3]) based on A mo( j i 
Ap^pfina^o, and Ap = pg na i. However, to estimate poste- 
riors for these means from the counted numbers, we must 
account for incompleteness and Poisson uncertainty. 

The prediction by S12 was for an ideal Kepler sam- 
ple complete out to 2 years (Subo Dong, private com- 
munication, 2012). To derive the expected number of 
super-eccentric proto-hot Jupiters in a sample of a lim- 
ited timespan P max we must account for incompleteness. 
Given that we can only detect proto-hot Jupiters with 
two or more transits, the completeness (with respect to 
this effect alone) C comp ranges from 100% at orbital pe- 
riods < P m ax/2 to 0% at orbital periods of P max . We 
update Equation Q to account for incompleteness: 
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Fig. 1. — Giant planets detected by non-Kepler transit surveys. All are Jupiter-mass (M p > 0.25Mj up ) planets orbiting stars with 
4500 < T c ff < 6500 K, logg > 4. The dashed lines are tracks of constant angular momentum corresponding to Pfinai = 3, 5, 10 days. The 
blue stars are moderately-eccentric Jupiters (-/V mo( j q). The circles have e < 0.2. The open symbols (black open circles and blue open stars) 
have 3 < P < 10 days and are used to calculate Np— p nna i o- The blue, diagonal striped region is where the predicted super-eccentric 
proto-hot Jupiters with e > 0.9 and 3 < P nna i < 10 days would reside. 



where the completeness C com p(e) is 

(^comp (^) 

1 , e < e(P max /2) 
(Pmax/Pfinai)(l - e 2 f' 2 - 1 , e(P max /2) < e < e mi 



(8) 
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pie, for P n 



1 - (2Pfinal/Pi 

2 years and Pg B 



s2/3 



1/2 



For exam- 



1.04 (Equation [§ to 0.868 (Equation |7]). ^ 
ness applies to the jBatalha et ah (2012 



3_days, r drops from 
This complete- 
ly ) sample, cov- 
ering two years, assumes that all Jupiter-sized planets 
exhibiting two transits in the Q1-Q8 data are detected, 
and neglects the effects of missing data. We revisit these 
assumptions in Section 4.3 and find that relaxing them 
does not account for the lack of super-eccentric proto-hot 
Jupiters. 

Next we describe the selection cuts we make to count 
-ZVmod.o (blue stars, Figure[T]), N P=Pfina ,ifl (open symbols, 
Figure [I]), and Np—psnai- Because the stellar parameters 
from the Kepler Input Catalog (KIC) are not reliable 
for stars outside the temp erature range 4500 < T c g < 
6500 K (Brown et al. 112011 ) we only include stars within 



this temperature range m both the Kepler and calibra- 
tion samples. We impose an aggressive cut of stellar sur- 
face gravity log g > 4 to exclude giant stars, because their 
KIC parameters are unreliable. From the Kepler sample, 
we remove three kn own false positives (KOI-425. 01,KOI- 
552.0 1, KOI-609.01, |Santerne et al.|2012[ see also |Morton| 
2012 1. Since it is not possible to detect planets along ex- 
actly the same angular momentum track, we follow S12 
and consider two Pfi na i intervals: 3 < Pfinai < 5 (Interval 
1) and 5 < Pg na i < 10 (Interval 2). We tabulate the 
counted numbers and their sources in Table [2] 

Each number of counted planets (Table [2]) is drawn 
from a Poisson distribution with an unknown mean. We 
wish to compute the expected number of super-eccentric 
proto-hot Jupiters using not the counted numbers but 
estimated posteriors for the mean numbers, incorporat- 
ing uncertainty. The probability of observing N p \ from a 
Poisson distribution with mean N p i is: 



prob(iVpi|7V pl ) 



exp[ 



-N 



pU 



(9) 



We estimate the posteriors of mean numbers (effectively 
inverting the Poisson distribution) via a Markov Chain 
Monte Carlo (MCMC) algorithm with the log-likelihood: 



JV p iln(iV p i)-hi(iV p i!)-iV : 



pi 



(10) 



and obtain iVpi posteriors, for which we include the me- 
dian and 68.3% confidence intervals in Table [2j These 
posteriors have medians slightly larger than the counted 
numbers because of the skewed shape of a Poisson dis- 
tribution at small values of the mean (JV p i < 10). It is 
more probable that we are observing fewer planets than 
the true mean number than vice versa. 

Next we compute the mean number of super-eccentric 
planets, N sup . To do so we insert the posteriors from 
Tablep] into Equation |3]) , making use of Equations ([6| 
and (TnT We perform this procedure separately for Inter- 
val 1 and Interval 2 and use P max = 2 year, obtainingn 
a iV sup posterior for each interval, which we combine to 



compute a total -/V sup (Figure 2k. The total expected 



' sup 

number is iV sup = 7^3 . This posterior represents a distri- 



bution of Poisson means. We transform the distribution 
of means into a distribution of expected values by sam- 
pling iV sup from iV sup according to Equation Each 
sample requires first drawinga mean (N sup ) from the 
distribution of means (Figure pi top row) and then draw- 
ing an observed number N sup from the Poisson distribu- 
tion with that mean. The expected observed number is 
7^3, with only a 0.5% probability of observing 



sup 

0. Therefore, despite small number statistics, we expect 
a readily detectable number of super-eccentric proto-hot 
Jupiters given the S12 assumptions. 

3.2. Prediction for transit light curve observables 

We expect to be able to identify super-eccentric proto- 
hot Jupiters in the Kepler sample by fitting their tran- 
sit light curves and identifying those for which the light 
curve model parameters are inconsistent with a circular 
orbit. A planet's orbital eccentricit y affects its transit 
light c urve in a number of ways (e.g. |Barnes|[2007] |Ford| 



et al. 2008 



Kipping 



20081). For long-period, highly ec 
centric, J uplter-sized planets, the most detectable effect 
is on the transit duration. For a wide range of periapse 
orientations relative to our line of sight, a planet on a 
highly eccentric orbit transits its star moving at a much 
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TABLE 2 

Counted planets 



e Interval Counted Mean a Sample b 

[days] 



0.2 < e < 0.6 


1 


3-5 


-Nmod,0 = 


= 5 


N mod,0= 6_ 2 


Cal 




2 


5-10 




= 3 


^ r mod,0= 3.7lj;g 


Cal 


unspecified 


1 


3-5 


A r p=Pfinal,0 = 


41 


A r P = PflnaL0= 45lg 


Cal 




1 


3-5 


Np = Pfinal = 


31 


A , P = Pfinal = 32l| 


Kep 




2 


5-10 


-?Vp=Pfinal,0 = 


= 9 


^P = Pfinal,0= 10 ± 3 


Cal 




2 


5-10 


-?Vp=Pfinal = 


19 


A r P = Pflnal= 20±f 


Kep 



a Median, with 68.3% confidence interval, of posterior of Poisson means, each 
defining a Poisson distribution from which the counted number may be sampled. 
b Kep = Kepler ; Cal = calibration non-Kepler transit-detected (Figure 111. 



Predicted mean 




Predicted number observed; 



10 15 
N supereccentric 



25 



Fig. 2. — Top: Predicted mean number of super-eccentric 
Jupiters (Interval 1: red dotted, Interval 2: blue dashed, total: 
black solid). Bottom: Sampling from above distribution of Poisson 
means to create a distribution of expected number observed. 

larger speed than if it were on a circular orbit with the 
same orbital period. For Jupiter-sized planets, one can 
distinguish the effects of the transit speed on the ingress, 
egress, and full transit duration from the effects of the 
transit impact parameter and/or limb-darkening, even 
with long-cadence Kepler data (DJ12). 



For each planet, we fit a Mandel & Agol (2002) tran- 
sit light curve model with the following parameters: the 
planetary-to-stellar radius ratio R p /R+, the orbital pe- 
riod P, the inclination i, the scaled semi-major axis 
a/i?*, and the quadratic limb darkening parameters pi 
and fi2- Assuming a circular orbit and a planetary mass 
much less than the stellar mass (M p << M*), one can 
recast the scaled semi-major axis a/i?* (by substitut- 

r 2i 1/3 

ing a = GM± (P/[2tt}) ) as the bulk stellar density, 

p+ = M±/(^ttR'1), which we will refer to simply as the 
stellar density hereafter. The stellar density measured 
from the light curve under the assumption of a circular 
orbit, Pcirc, is related to the true stellar density, p* by: 



P*{e,u)g (e,w) = p circ , 

where 
1 + e sin w 

g{e,u}) = 



vT 



(11) 



(12) 



is approximately the ratio of the observed transit speed 
to the transit speed that the planet would have if it were 



on a circul ar orbit with the same orbital period (see Kip- 
ping|2010 and DJ12 for a detailed derivation). The argu- 



ment of periapse u represents the angle on the sky plane 
(lo = 90° for a planet transiting at periapse). 

We determine p c i re by fixing e = 0, allowing the stel- 
lar density to vary as a free parameter in the light curve 
model. The resulting p C i rc is determined entirely by the 
shape and timing of the light curve. We then compare 
p c i rc to the value of p* determine through other meth- 
ods (i.e. stellar models fit to the temperature and sur- 
face gravity determined through colors or spectroscopy). 
Although g is degenerate with the host star's density 
(Equation 111, a loose (order-of-magnitude) constraint 
on p* is sufficient for a tight constraint on the eccen- 
tricity (DJ12), measurement of which we will describe 
and perform in Section [4] For now, we use p CKC /p*- If 
Pcirc/p* is very large, then g must be large, and there- 
fore the planet is moving more quickly during transits 
than a planet with orbital period P on a circular orbit. 
In Appendix [Bj we summarize how our approach avoids 
problems caused by incorrect stellar parameters. 

3.2.1. Expectations for super- eccentric planets 

We perform a Monte Carlo simulation to predict the 
signature in the transit light curve observable Pctrc/p* ex- 
pect ed f rom the super-eccentric proto-hot Jupiters (Sec- 
We generate two-dimensional (2D) probability 



3.1) 



tion 

distributions in (P, p C h C /p±) in Figure [3j where P is the 
orbital period, as follows: 

1. We begin with an assumed Pfinai- 

2. Using the completeness Equation |8j we generate a 
distribution of eccentricities {e{\ with a normaliza- 
tion constant C n orm following: 



Prob(e) 







0.9 < e 



or e < 

< ' infix 



0.9 



(13) 



3. For each eccentricity, we compute the correspond- 
ing orbital period Pi and randomly select an argu- 
ment of periapse uji. Assuming a Sun- like star, we 
compute the scaled semi-major axis a^/i?^. 



4. We compute the transit probability: 
R± 1 + e, ; sin 



prob 



transit 



1 



(14) 
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Then we select a uniform random number between 
and 1. If the number is less than the transit 
probability, we retain (ej,u>i) in the distribution. 



5. Then we compute Pdic/p* using Equation (11) 



We use the above procedure to generate four plots, cor- 
responding to different Pfi na i (Figure [3]) . In the fourth 
panel, instead of using a single Pfinai, we draw the P nna i 
of each trial from the observed Ap = p nna i, weighting each 
Pfinai by iV mo d,o/Np=Pfinai,o in the two intervals. We 
see that a population of super-eccentric Jupiters will 
manifest itself as a collection of light curves with as- 
trophysically implausible p C j rc of 10-1000 times the es- 
timated values for p+. The super-eccentric proto-hot 
Jupiters will have orbital periods that range from P = 
3 days/(l — 0.9 2 ) 3 / 2 = 36 days to the observational times- 
pan of two years. About 90% of the expected planets 
have Pcirc/p* > 10, making them easy to identify. 

3.2.2. Proto-hot Jupiters with 0.6 < e < 0.9 

S12 focused their prediction on super-eccentric planets 
with e > 0.9. However, we also expect to find proto-hot 
Jupiters with less extreme eccentricities (0.6 < e < 0.9) 
along the same Pg n al track. We repeat the procedure 
in I33TT1 for the interval 0.6 < e < 0.9. The over- 
all occurrence rate for this interval is 0.61 relative to 
■^mod- As shown in Figure [ij the proto-hot Jupiters in 
the 0.6 < e < 0.9 range have shorter orbital periods 
(6 < P < 121 days). However, their transit durations 
and the inferred stellar density from a circular fit are not 
as strikingly anomalous as for the super-eccentric proto- 
hot Jupiters, making them less easy to identify. There- 
fore, we do not focus on these objects but discuss them 
further in the conclusion (Section 6). 

4. RESULTS: A PAUCITY OF PROTO-HOT JUPITERS 

We search for the super-eccentric proto-hot Jupiters 
predicted by S12 and find significantly fewer than ex- 
pected. We describe our s earch procedure and present 
our measurements (Section 4.1 |, assess the significance 
of this null result (Section 3T2"), and consider wh ethe r 
observational causes are likely responsible (Section 4.3). 



4.1. Transit light curve observables for potential 
proto-hot Jupiters 

We begin by identifying planet candidates that con- 
form to our selection criteria using the Kepler Planet 
Candidate Data Explo rer 6 , an interactive t able of the 

with a de- 
Applying 



candidates reported by Batal ha et al.|(|2 012), 
sign similar to the EOD ( | Wright et al.pOTT| 
the same criteria as in Section 2, we identify candidates 
with R > 8P® and stellar parameters 4500 < T c g < 
6500 K, log g > 4. We restrict the orbital periods to 
those between 36 days and 2 years, with the upper limit 
chosen to match S12 for ease of comparison and the lower 
limit corresponding to Pfi na i = 3 for e = 0.9. We are left 
with 41 planet candidates, including KOI-1474.01 (D12). 

For each candidate, we retrieve the Q0-Q9 publicly- 
available dat a from MAST. We extract the transits us- 
ing AutoKep (|Gazak et al.|2012[ ) and perfor m an MCMC 
fit using the l'ransit Analysis Package (TAP;|Gazak et alT 



2012|. We fix e — but allow all othe r parameters to 
vary, i ncluding noise parameters for the |Carter fc Winn| 
( 2009 \ wavelet likelihood function and first-order polyno- 
mial correction terms. We use short-cadence data when 
available. We obtain each candidate's p c i rc posterior. 
Next we follow 7 Section 3.3 of D 12 to compute a p+ 



posterior for each host star using the Takeda et al. ( 2007 1 



stellar evolution models and the estimated effective tern 
perature an d surface gravity liste d in the Kepler candi- 
dates Table ( |Batalha et al.|[2~012[ ). For stars with stellar 
properties measured from spectroscopy, we adopt con- 
servative uncertainties of 100 K for T e g and 0.1 for log q 



(all a bove the uncertainties reported in Buchhave et al 
2012 1. For stars for which the effective temperature and 



surface gravity were derived from the KIC colors, we 
adopt uncertainties of 200 K for the effective temper- 



ature and 0.4 for log<? ( Brown et aLpOll |. We describe 
exceptions to this procedure, as well as additional cuts 
that left us with 31 planet candidates, in Appendix |C| 

Finally, we combine the p c - uc and posteriors into 
a posterior of Pdic/p* for each candidate, marginalized 
over all other parameters. In Figure[5l we plot the result- 
ing values on top of the probability distribution for pre- 
dicted super-eccentric proto-hot Jupiters (Figure |3j panel 
4). None of the candidates fall in the high-probability 
area of the prediction. We indicate candidates with 
known companions in their system with blue bars; none 
can have e > 0.9 and 3 < P nn ai < 10 days without its or- 
bit crossing a companion's. As expected, all 8 candidates 
with companions have p C irc/p* close to 1. 

Two candidates without companions, KOI-1474.01 
(D12) and KOI-211.01 both have p c -„c/p* > 10. 
However, KOI-211.01 has a long orbital period of 
372.108 days, requiring e > 0.954 to have P nna i < 
10 days, corresponding to p c i xc / P* > 277 for a transit 
at periapse. This requirement is reflected in the slope of 
the high-probability density region shown in Figurep](i.e. 
the longer the orbital period, the larger the predicted 
Pcirc/p*)- The probability of KOI-211.01 having e > 0.9 
and 3 < P nna i < 10 days is 12%. D12 found that KOI- 
1474.01 has e = 0.81+g ^ and P final = U±% days. The 
probability of it having e > 0.9 and 3 < P nna i < 10 days 
is also 12%. In assessing the consistency of the observa- 
tions with the prediction of Section |3j we will fully con- 
sider the possibility that KOI-211.01 and KOI-1474.01 
might be super-eccentric proto-hot Jupiters. 

We emphasize that it is not surprising that many of 
the candidates lie in the low-probability region (likely 
Jupiters with e < 0.9 or Pfi na i > 10 days, of which there 
may be any number). It is only surprising that we do 
not see half a dozen in the high-probability region. 



4.2. Statistical significance of lack of proto-hot Jupiters 

None of the observed candidates lie in the high- 
probability-density region of Figure [5] so i t is unlikely 
that the half a dozen predicted (Section 3.1) super- 



planetquest . jpl .nasa.gov/kepler 



7 Instead of imposing a prior on the stellar mass, metallicity, and 
age from a TRILEG AL (TRIdimensional modeL of thE GALaxy; 
|Girardi et al.|[2005| synthetic Kepler field population, we assume 
a uniform prior on these model paramet e rs, be cause a similar prior 
was already imposed by |Batalha et al.| (|2010[ l to generate the ef- 
fective temperature and surface gravity in the KIC. 

8 KOI-433.02 has a low p c irc/p*- We discuss it further in Ap- 
pendix [C] 
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Fig. 3. — Top: 2D posterior, orbital period P vs. p c irc/p*> for planets with e > 0.9 and P nIla i = 3, 5, 10 days (panels 1-3) or P nIla i drawn 
from Kepler hot Jupiters with 3 < P < 10 days (panel 4). Bottom: Posterior p c i rc /p* marginalized over orbital period. Proto-hot Jupitcrs 
with e > 0.9 should have anomalously large p c i rc measured from the transit light curve compared to their estimated p* , making them easy 
to identify. We expect half a dozen super-eccentric proto-hot Jupitcrs in the high probability density region. 
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Fig. 4. — Top: 2D posterior for orbital period P vs. p c irc/p* for planets with 0.6 < e < 0.9 and Pfl na i = 3,5,10 days (panels 1 - 3) or 
with Pfi Ila i drawn from Kepler hot Jupiters in the interval 3 < P < 10 days (panel 4). Bottom: Posterior p c irc/p* marginalized over orbital 
period. Proto-hot Jupiters with 0.6 < e < 0.9 do not typically have such large p c i rc /p* as their super-eccentric (e > 0.9) counterparts 
(Figure making them less easy to identify. 



eccentric proto-hot- Jupiters are present but missed. If we 
were certain that none of the candidates has e > 0.9 and 
3 < -Pfinal < 10 days, the probability that observed num- 
ber of super-eccentric proto-hot- Jupiters agrees with the 
prediction would simply be 0.5%. (This is the probabil- 
ity, computed in Section 3, of observing super-eccentric 
proto-hot Jupiters given the Poisson uncertainties in the 
observed number of super-eccentric Jupiters and in num- 
bers used to compute the prediction.) However, there is a 
small chance that there are indeed super-eccentric proto- 
hot-Jupiters among the sample but that they just so hap- 
pen to have their periapses oriented in the narrow range 
of angles producing an unremarkable Pdrc/p*- Therefore 
we use a Monte Carlo procedure to assess the consistency 



of Pcirc/p* posterior derived for each candidate with the 
predicted population of super-eccentric planets. 

We first use the p C uc/ P* posteriors to generate an ec- 
centricity posterior for each candidate, via a MCMC ex- 
ploration of a limited set of parameters: p c \ rc , p±, e and u> 
(as outlined in DJ12, Section 3.4). Although we can only 
make a tight eccentricity measurement when the planet's 
eccentricity is large (DJ12), the broad eccentricity poste- 
rior for the typical candidate here is useful for this pur- 
pose: it contains very little probability at the high eccen- 
tricities corresponding to e > 0.9, 3 < Pfinal < 10 days. 
We then perform 10 6 trials in which we randomly select 
an eccentricity from each candidate's eccentricity poste- 
rior. We compute Pfi na i and count N sup in Intervals 1 
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Fig. 5. — Expected 2D posterior for orbital period P vs. p c irc/p* 
(taken from panel 4 of Figure py . The values we measured for 
our 31 candidates are overplottea. Thin, blue bars: candidates 
with companions in their systems. Thick, red bars: candidates 
with no known companions. We do not see the expected half a 
dozen candidates in the region of high-probability density. Blue 
diamonds: KOI-1474.01 (bottom), KOI-211.01 (top). 

and 2. If both are greater than or equal to the respec- 
tive numbers drawn from posteriors in Figure [2j bottom 
panel (red dotted and blue dashed curves), we count the 
trial as a success, meaning that at least as many super- 
eccentric Jupiters as predicted were detected. 98. 7% 9 of 
trials were unsuccessful. We exclude the candidates with 
known companions from this procedure (Figure [5] thin 
blue bars) , because it so happens that none of them can 
have e > 0.9 and 3 < P nn ai < 10 days without cross- 
ing the orbit of another candidate in the system. We 
find that, with 98.7% confidence, we detected too few 
super-eccentric proto-hot Jupiters to be consistent with 
the prediction of Section [3] For example, 53% of trials 
had super-eccentric proto-hot Jupiters, 87% had 1 or 
fewer, and 98% had 2 or fewer. From these trials, we 



measure a N sup posterior with a median N s , 



4.3. Observational explanations for a paucity of 
proto-hot- Jupiters 

As discussed in Section [2~2} observational assumptions 
affect the predicted number of super-eccentric proto-hot 
Jupiters, which could be could be too large if the sample 
of Kepler hot Jupiters were to contain a large fraction 
of false positives or if the sample of long-period Kepler 
Jupiters were incomplete due to missing data or noise. 
Repeating the calculations in Section |3.1| we find that 
55% of our Kepler hot- Jupiter sample - from which we 
already removed known false-positives - would have to 

9 We note that if we treat Intervals 1 and 2 combined instead 
of separately, the percentage of unsuccessful trials decreases some- 
what to 97.9%. However, as the prediction applies to both intervals 
separat ely, w e argue that the separate treatment is appropriate. In 
Section |5.1| we consider the possibility that only hot Jupiters in 
Interval 1 have undergone tidal circularization. 



D. 
Q. 
Q 

O 



1000 r 



100 r 



10 r 




10 

Period (days) 

Fig. 6. — Combined differential photometri c pre cision (CDPP) 
for Jupiter-hosting stars using cuts from Section [3.1| A large CDPP 
indicates more noise. Triangles: Hot Jupiters used to compute 
iVp = pfl na i. Circles: Sample of potential super-eccentric Jupiters. 
Diamonds: Other Jupiters. The dashed line, CDPP = 300, marks 
a conservative completeness cut-off for the CDPP out to 2 years. 

be false positives in order for our confidence that we 
observed fewer than expected super-eccentric proto-hot 
Jupiters to drop from 98.7% to 90%. Such a high false- 
positive rate is unlikely given the low a priori likelihood o f 
false positives ( |Morton k. Johnson|2011b| |Morton|2012[ ). 

In Section [3j we accounted for incompleteness due to 
the two year timespan of the data. We neglected the 
effect of missing data. In Appendix |Dj we describe our 
procedure for estimating the incompleteness due to miss- 
ing quarters. We repeat the procedure of Section 4J2 
and find that our confidence only d ecreases slightly, from 
98.7% to 98.5%. Additionally, the |Batalha et al.f ( |20T2 1 
Kepler candidate collection may be less complete than 
we assumed for planets whose second transits occur in 
Q7 or Q8. The Kepler pipeline was run on the Q1-Q6 
data, and then the planetary parameters were fit using 
Q0-Q8 data. Candidates exhibiting only one transit in 
the Q1-Q6 data were identified, but possibly with less 
completeness. We repeat the procedure of Section |4.2| 
with P max = 1.5 years instead of two years, imposing 
an overly conservative requirement that a candidate ex- 
hibit two transits in Q1-Q6. We also incorporate miss- 
ing data (Section [D| . Our confidence only decreases to 
97.2%, and therefore these effects do not account for the 
lack of super-eccentric proto-hot Jupiters. 

Although the planets we considered here are relatively 
easy to detect due to their large transit depths, possi- 
bly the sample is incomplete around the noisiest stars. 
In Figure [6j we plot the combined differential photomet- 
ric precision (CDPP) for all stars with 4500 < T cS < 
6500 K, log g > 4 hosting Jupiter-sized candidates. A 
higher CDPP indicates a noisier star. We consider how 
the expected number of proto-hot- Jupiters would decline 
if we assumed the sample were only complete for CDPP 
< 300 (dashed line). If we restrict the sample to Jupiters 
with CDPP < 300, we are left with 23 Kepler Jupiters 
with 3 < P < 5 days and 12 with 5 < P < 10 days. 
We repeat the calculations of Section [3] and find that the 
expected number of super-eccentric proto-hot Jupiters 
would be 5+3. The chance of observin g w ould be 1.8%. 
We repeat the procedure of Section 4J2 and find that 
we nonetheless detected fewer super-eccentric proto-hot 
Jupiters than expected with 96% confidence. 

Therefore, false positives and incompleteness are un- 
likely to account for the paucity of proto-hot Jupiters in 
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the Kepler sample. However, we recommend revisiting 
our calculations as new false positive estimates and new 
candidates become available. 

5. EXPLAINING THE PAUCITY OF PROTO-HOT 
JUPITERS 

So far (Sections 1-4) we have been considering a sce- 
nario in which hot Jupiters begin beyond the ice line 
on super-eccentric orbits — caused by gravitational per- 
turbations from a companion (e.g. stellar binary Kozai, 
planetary Kozai, planet-planet scattering, secular chaos) 
— and subsequently undergo tidal circularization along 
a constant angular momentum track, reaching a final 
orbital period Panai- This process is known as high- 
eccentricity migration (HEM). We schematically summa- 
rize this (black arrows) and other possible origins for hot 
Jupiters (white and gray arrows), as well as moderately- 
eccentric Jupiters with 3 < Pfi na i < 10 days, in Fig- 
ure [7] The correspon ding populations fro m the RV- 



Wright et al. 20121 are plotted 



41). In Sections |5J 
that Jupiters began 



detected sample (EOD, 

in Figure [8] Now we relax previous assumptions about 
HEM (Section 2.2) and explore how we can account for 
the lack of supe r-ec centr ic p roto-hot Jupiters (Section 
and 5.2 we relax the assumption 
eyond the ice line, finding that this 
possibility could indeed account for the lack of super- 
eccentric Jupiters. In Section |5.3| we consider the par- 
ticular case of HEM via the Ko zai m echanism with a 
planetary perturber. In Section |5.4[ we relax the as- 
sumption of a steady current of hot Jupiters produced 
by HEM but find that a lack of steady current is un- 
likely to account for the lack of super-eccentric proto-hot 
Jupiters. In Section [B~5| we place an upper-limit on the 
fraction of hot Jupiters caused by Kozai perturbations 
from a stellar binary companion. 

5.1. No tidal circularization: hot Jupiters and 
moderately- eccentric Jupiters implanted interior to 
the ice line 

Rather than starting on highly-eccentric orbits exte- 
rior to the ice line, hot Jupiters and moderately-eccentric 
Jupiters may have been placed directly into the region we 
observe today. The moderately-eccentric Jupiters (blue 
stars, Figure [I] [7| [8| observed along the angular mo- 
mentum tracks have may not have undergone tidal cir- 
cularization but may have been placed there by whatever 
mechanism eccentric Jupiters interior to the ice line (gray 
region, Figures [7] and fsj) . This underlying population 
could originate from planet-planet scat tering or secular 
chaos, or possibly from disk migration. |Wu fc Lithwick 
( |2011[ ) found that secular chaos should produce a num 
ber of moderate-to-high ecc entricity "warm Jupiter s" in 
the region from 1 to 1 AU . |Goldreich fc Sari| ( |2003| ) and 
Sari & Goldreich (2004) argued that disk migration can 
potentially excite moderate eccentricities through reso- 
nance torques, but recently Dun hill et al.|[2012| mod- 
eled planet-disk interactions using high-resolution three- 
dimensional simulations and found that disks are unlikely 
to excite the eccentricities of giant planets. 

In Figure [5J the blue stars look as if they could be an 
extension of the distribution in the gray region. In Ta- 
ble [3j we compute the occurrence rate of giant planets in 
the RV-discovered sample in different regions of Figure[8j 
The occurrence rate per log semi-major axis interval of 



TABLE 3 

OCCURRENCE a OF JUPITERS DETECTED BY RV SURVEYS 



Period range 
(days) 


Eccentricity 


Count 


Poisson 
range b 


Number 
per log 10 a 


3 < Pflnal < 10 

10 < P < 250 


0.2 < e < 0.6 
0.2 < e < 0.6 


2 

18 


r, 7 +2.0 

z - ' -1.3 

19±| 


8t 6 4 
20t| 


3 < P < 5 
3 < P < 10 
5 < P < 10 
5 < P < 250 
10 < P < 250 


e < 0.2 
e < 0.2 
e < 0.2 
e < 0.2 
e < 0.2 


13 
17 

1 

27 
23 


14±i 
18±! 

4 7+2-5 
*■ '-1.8 

28+ 6 5 
24 ±5 


23+1 2 

25+J 

25t 6 5 



a Numbers do not account for RV observational biases 
b Range of Poisson means from whose distributions the count could 
have been drawn. The median is higher than the counted number 
due to the skewed Poisson distribution shape. 



moderately-eccentric Jupiters with 3 < Pg na i < 10 days 
is less than or equal to that in the gray region (10 < P < 
250 days). Therefore a separate mechanism for producing 
the blue stars apart from direct implantation may not 
be necessary. If non-tidal implantation was dominant, 
the number of moderately-eccentric Jupiters should not 
be used to predict the number of super-eccentric Jupiters 
because the moderately-eccentric Jupiters did not tidally 
circularize from super-eccentric orbits. 

If moderately-eccentric Jupiters with 3 < Pfinai < 
10 days did not undergo tidal circularization, hot 
Jupiters themselves could be part of a continuous dis- 
tribution of circular Jupiters interior to the ice line (Fig- 
ure lf\ and pi black region) , which must have migrated 
somehow. Disk migration effectively produces planets 
on circular orbits, but seems inconsistent with the high 



obliquities of hot Jupiters orbiting hot stars (Winn et 
[aL] 120101 |Albrecht et al.||2012[ but see also |Rogers et 
^12f. ' 



|ar p012| ). However, disk migration may have produced 
some or all of the well-aligned hot Jupiters, if their low 
obliquities are not the result of tidal realignment. 

If the cut-off for tidal circularization is 3 days, rather 
than 10 days, Jupiters on circular orbits with P > 3 days 
would actually be part of the so-called "period-valley," 
rather than the hot Jupiter pile-up. The period valley 
refers to the region exterior to hot Jupiters but interior 
to the ice line (P < 250 days), where giant planets are 
scarce. The divide between hot Jupiters and the period 
valley (i.e. if it is 3 days, 10 days or some other va l ue) is 
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et al.||2003| |Wright et aL]|2009| Wittenmyer et al.j [2UT 
The observed "edge" of hot Jupiters in ground- bas 



ft 



upiters in ground- based 
transit surveys may be partially caused by a combina- 
tion of the reduced geometric transit probability of long- 
period planets and inefficie ncy of ground-based transit 



surveys in detecting them (Gaudi et al. 2005), rather 
than a drop in the intrinsic occurrence rate. We note 
that the distribution of giant planets infe rred from the 
Kepler Mission, as sessed out to 50 days ( Youdin 2011| 
Howard et al.|2012 ) , has no such edge. However, Wright 



et al. ~ p009 [) detect an edge at approximately 0.07 AU 
•s) in 



(5 days) in their RV survey, which suffers from different 
(but less severe) biases than transit surveys. Thus, the 
existence and location of the cut-off remains uncertain. 
Therefore the cut-off may in fact be between P = 
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Fig. 7. — Schematic of pathways (arrows) for creating the observed giant planet populations, which we assume formed beyond the ice 
line and reached semi-major axes interior to the ice line via the stellar binary Kozai mechanism, planetary Kozai mechanism, planet-planet 
scattering, secular chaos, or disk migration. "Track" refers to the angular momentum range under consideration, i.e. 3 < P nna i < 10 days. 
The black arrows indicate the path that we have assumed throughout the paper for HEM caused either by a planetary or stellar perturbcr. 
For example, a Jupiter may be perturbed by a stellar binary companion, follow the black arrow to the region of super-eccentric Jupiters with 
3 < Pfinai < 10 days (blue, diagonal stripe region), undergo tidal circularization along its angular momentum track to e < 0.9 (horizontal 
orange striped region), become a moderately-eccentric Jupiter (blue stars), and eventually achieve hot- Jupiter- hood (red triangles). The 
purple dashed arrows represent planetary Kozai, which (depending on the perturber) may cause the proto-hot Jupiter to undergo low- 
eccentricity excursions off the track (Smadar Naoz, private communication, 2012). The other color arrows indicate alternative pathways 
caused by secular chaos or scattering (white) or disk migration (gray), and colors and patterns of the boxes correspond to the regions of 
parameter space in Figure [8] See text for detailed discussion of the possible scenarios represented here. 
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Fig. 8. — RV-discovered planets from the EOD | |Wright et al.|2012| queried on September 6th, 2012) with M p sini > 0.25M Ju p and stellar 
properties 4500 < T e ff < 6500 K, logg > 4. Note tnat tne planets from Figure [l] are not included here because they were discovered by 
the transit method. The dashed lines represent tracks of P nna i = 3,5, 10 days. The shaded and patterned regions correspond to Figure n\ 
Within the 3 < P nna i < 10 days angular momentum tracks are hot Jupiters (red triangles), moderately-eccentric Jupiters with 0.2 < e < OS 
(blue stars), Jupiters with 0.6 < e < 0.9 (orange horizontal striped region), and super-eccentric Jupiters (blue, diagional-striped region), 
The RV-discovered planet HD 17156 b lies in the orange, horizontal striped region, and the RV-discovered planet HD 80606 b lies in the 
blue, diagonal-striped region. Period valley: grey region denotes Jupiters with P nna i > 10 days but interior to the ice line, and the black 
region houses circular Jupiters interior to the ice line but exterior to hot Jupiters. 
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3 days and P = 10 days. In Figure [S] the edge of the pile- 
up of circul ar Jupiters appears t o end at around 0.057 AU 



(5 days), as Wright et al. (2009) found. If we separate the 



hot Jupiters below this cut-off, we recover a pile-up of hot 
Jupiters: in the region from 3 < P < 5 days, we observe 
an excess of circular Jupiters inconsistent with the occur- 
rence rate in the period valley by a factor of 3 (Table [3]) . 
If the cut-off for hot-Jupiters is truly 5 days, the pre- 
diction for super-eccentric proto-hot Jupiters should be 
based only on the number moderately-eccentric Jupiters 
with 3 < Pfinai < 5 da ys. In that case (repeating 



the calculations of Section 3.1), we expect to find only 



3^2 super-eccentric proto-hot Jupiters with e > 0.9 and 
3 < -Pfinai < 5 days, and our con fidence that we found 
fewer than predicted (Section 4.2) drops to 86%. 



5.2. Some or all proto-hot Jupiters may have bypassed 
the e > 0.9 portion of the P nna i track 

Alternatively, the typical hot Jupiter may have under- 
gone tidal circularization but bypassed the high eccen- 
tricity phase, starting on the HEM track with 0.6 < 
e < 0.9 in the region indicated by orange stripes in 
Figure [7] and |8] (or even in the 0.2 < e < 0.6 region). 
For Pfinai < lTJ days, a Jupiter would begin the HEM 
track at an orbital period less than 120 days, or 0.5 AU 
around a Sun-like star. The Jupiter is unlikely to have 
formed here — the critical core mass required to accrete 
a massive atmosphere most likel y exceeds the amount 
of refractory materials available (Rafikov 2006) — but 
may have been delivered to this region via planet-planet 
scattering or secular chaos. Assuming a steady-flux of 
proto-hot Jupiters into the orange striped region, if the 
two moderately-eccentric, RV-detected 10 Jupiters (cor- 
2.7±f° 3 ) in Figure Is] (blue stars) 



responding to N mo( 

originated from the orange region, we woulcl expect 
A^mod / C , comp (e)|e|" 1 de I / / C com p(e)|e| _:L de 



0.6 



0.2 
2.0 



2.7Zt% x °- 616 = !-7 



+1.2 
-0.8 



proto-hot Jupiters in the orange region. We indeed see 
one such planet, HD 17156b (Figure [8]). 

Since we observe HD 80606 b in the blue striped re- 
gion, all proto-hot Jupiters would not necessarily begin 
in the orange striped region. Planet-planet scattering or 
secular chaos may place proto-hot Jupiters in both the 
orange striped region and the blue striped region, with 
the majority in orange striped region. The proto-hot 
Jupiters in these two regions would be created by the 
same dynamical processes responsible for Jupiters with 
-Pfinai > 10 days (failed hot Jupiters), of which we observe 
more with 0.6 < e < 0.9 than with e > 0.9 (though this 
may be partly due to observational bias). The overall 
picture of this scenario is that proto-hot Jupiters start 
the HEM track interior to the ice line with eccentricities 
similar to those of planets we observe in the period valley, 
rather than starting with e — > 1 beyond the ice line. 

10 In this calculation, we use the RV-detected sample. Even 
though transit probability is constant along an a nna i track, ground- 
based transit survey are still strongly biased against detecting plan- 
ets transiting with longer orbital periods. RV samples suffer from 
their own biases against long period and eccentric planets, which 
we do not account for here. 



5.3. Proto-hot Jupiters created by planetary Kozai 

Proto-hot Jupiters created by distant stellar compan- 
ions likely decoupled by time they reach P < 2 years. 
However, those created by nearby planetary companions 
may still be coupled to their perturbers and spend much 
of their evolution at low e (Smadar Naoz, private commu- 
nication, 2012). The possibility that a proto-hot Jupiter 
spends part of its time off the HEM track due to Kozai 
oscillations is indicated by white, purple dashed arrows 
in Figure [7] We clarify that the possibility that superec- 
centric proto-hot Jupiters spend time at low eccentrici- 
ties does not reduce their expected number. S12 demon- 
strate that tidal dissipation primarily occurs during high 
eccentricity intervals. Regardless of how much time the 
Jupiter spends off its a nna i track during low-eccentricity 
Kozai phases, it spends the same total amount of time 
on the afinai track undergoing tidal dissipation. Therefore 
Equation ([3]) predicts the total number of super-eccentric 
Jupiters observed on the track. 

However, there are two ways in which the plan- 
etary Kozai could cause fewer super-eccentric proto- 
hot Jupiter than computed in Section 3. First, even 
the moderately-eccentric calibration proto-hot Jupiters 
( A^mod.o blue stars in Figure [I]) could potentially still be 
coupled to a nearby planetary perturber (Smadar Naoz, 
private communication, 2012). We note that several of 
the moderately-eccentric Jupiters (HAT-P-34 b, HAT- 
P-31 b, and WASP-8 b) have lin ear trends in the RV 
observations ( |Wright et al. 2011), indicating the pres- 
ence of a companion in the system; it may be the com- 
panion is sufficiently massive, inclined, and nearby to 
remain coupled. If the moderately-eccentric Jupiters are 
still undergoing Kozai oscillations, we may be observing 
them in the low-eccentricity portions of their cycles, and 
they may actually be tidally dissipating on a track with 
-Pfinai < 3 days. If so, they should not have been used to 
compute the expected number of super-eccentric proto- 
hot Jupiters with 3 < Pfinai < 10 days. We clarify that 
even if the moderately eccentric Jupiters oscillate to high 
eccentricity intervals corresponding to a dissipation track 
with Pfinai < 3 days, they should not be used to compute 
the number of super-eccentric Jupiters with Pfinai < 3 
days, because Equation ^ only applies to planets ob- 
served on the track. 

Second, assuming that the proto-hot Jupiter and its 
planetary perturber formed co-planar in the disk, they 
likely underwent scattering to achieve the mutual incli- 
nations necessary for Kozai. This scattering process may 
have delivered the proto-hot- Jupiter interior to the ice 
line, leading to the scenario described in Section |5~2| in 
which the proto-hot Jupiter embarks on its HEM track at 
a < 0.5 AU and e < 0.9. Indeed, in planetary Kozai, the 
proto-hot Jupiter can sometimes only reach a maximum 
eccentricity of only 0.6 < e < 0.9, rather than e —> 1 
(Smadar Naoz, private communication, 2012). 

5.4. Alternatives to the "steady current" approximation 

The S12 prediction of a readily observable number 
of super-eccentric proto-hot Jupiters assumed a "steady 
current" of proto-hot Jupiter production. We only ex- 
pect a steady current if: a) the rate of hot Jupiter pro- 
duction throughout a star's lifetime is constant (e.g. that 
a hot Jupiter is just as likely to be produced between 4.1- 
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4.2 Gyr as it is during the first 100 Myr), or b) although 
the hot Jupiter production rate may, for example, be 
restricted to early in a star's lifetime, new planetary sys- 
tems are being produced at a constant rate throughout 
the Galaxy. In the latter case, we would see proto-hot 
Jupiters along all aa aa \ tracks only if the sample included 
stars in the stage of hot Jupiter production. Nonetheless, 
even if neither is true, we would still expect to observe 
proto-hot Jupiters. However, they would restricted to 
the narrow range of a nna i tracks for which the circular- 
ization timescale is of order a star's age, instead of being 
found along all a nn ai tracks in proportion to the circular- 
ization timescale (e/e cx a nnal , Equation pi) . 

Assumption (a) seems unlikely. In the HEM mech- 
anisms proposed (Section 1), proto-hot Jupiters are 
spawned on instability timescales (planet-planet scatter- 
ing, secular chaos) or the Kozai timescale, which are 
unlikely to always coincide with the typical stellar life- 
time. More likely, the distribution of timescales is uni- 
form (or normal) in order of magnitude and thus most 
proto-hot Jupiter s are spawned early in their host stars' 



lifetimes. Indeed, Quinn et al 



(20121 recently discovered 



hot Jupiters in the bOO Myr lieehive cluster and found 
that the, accounting for the cluster's enhanced metallic- 
ity, the hot Jupiter occurrence rate is consistent with that 
of the solar neighborhood. Regarding assumption (b), in 
practice we expect young stars to be rotating too rapidly 
to be amenable to Doppler observations and too uncom- 
mon in our stellar neighborhood to make up a represen- 
tative sample of transit surveys. Therefore the steady 
current approximation is unlikely to hold. 



In spired by popu lation simulations by Hansen (20101 
and [Hansen ( 2012[) , we simulate an extreme scenario m 
which every proto-hot Jupiter in the observable sample is 
created simultaneously (Figure [9) left panel). We be em 
with a population of Jupiters uniformly distributed in 
eccentricity and semi-major axis, extending to 10 AU; 
gray open circles had initial semi-major axes interior to 
1 AU (representing the possibility that planets can begin 
HEM interior to the ice line, as discussed in Section 5.2). 



Then we evolve the tidal evolution equations (Equations 

fand the corresponding a /a) until Jupiters with Pfi na i < 
days have circularized. We overplot tracks of constant 
angular momentum (dotted lines) , as well as lines defined 
by a constant "orbital change timescale," 



t, 



= (a/a) 2 + (e/e) 2 



-1/2 



(15) 



(dashed orange lines) which match the ag na i tracks at 
low eccentricities. Although there is no steady current, 
we see a "track" consisting of a) Jupiters along the same 
afinai track but with different starting eccentricities/semi- 
major axes, and b) Jupiters along close, adjacent a nna i 
tracks (those along the slightly larger a nna i track have 



'final 



higher eccentricities because e/e 

For comparison (Figure p)J right panel) we perform a 
simulation featuring a steady current of hot Jupiters. In 
this case, the proto-hot Jupiters are distributed over a 
range of angular momenta tracks but are most common 
(black diamonds) along the track where the tidal circu- 
larization time is order the total elapsed time (i.e. most 
of the proto-hot Jupiters began their HEM early in the 
lifetime of the oldest stars in the sample). The smaller 
ifinai tracks (red squares, blue X) are more sparsely pop- 



ulated because these planets circularize very quickly and 
we just happen to be catching some. The left and right 
panels are not strikingly different. Particularly with a 
small observational sample size, we are unlikely to be 
able to be able to distinguish whether we are seeing a 
narrow range of a nna i tracks due to a lack of steady cur- 
rent (left) or simply due to a higher relative population 
along the a nna i track of order a stellar age (right). 

However, without a steady current of hot Jupiters, the 
number of moderately-eccentric Jupiters in a Pfi na i range 
cannot be used in Equation ^ to predict the number of 
super-eccentric Jupiters. The number of super-eccentric 
Jupiters would depend on the initial conditions generated 
by HEM mechanisms: the relative number of hot Jupiters 
along adjacent angular momentum tracks and beginning 
at different semi-major axes along the same track. How- 
ever, if the initial eccentricities are roughly independent 
of semi-major axis, the distribution today would resem- 
ble one resulting from a steady current. The observed 
lack of super-eccentric proto-hot Jupiters would require 
very fine-tuned initial conditions, such as substantially 
fewer Jupiters beginning along a slightly larger angular 
momentum track. Therefore the paucity we found in Sec- 
tion 4 is unlikely to be fully accounted for by a lack of 
steady current. 

Throughout this work, we have assumed that proto-hot 
Jupiters can travel along HEM tracks with circulariza- 
tion timescales of order a stellar lifetime. However, a 
nearby planetary perturber can potentially permanently 
remove the proto-hot Jupiter from the angular momen- 
tum track before tidal circularization decouples it. For 
example, in secular chaos, a planet may be perturbed to a 
high eccentricity, begin to circularize along a track with 
a timescale longer than the chaos timescale, but then 
be chaotically perturbed by nearby planets again. If all 
proto-hot Jupiters are created early in the star's lifetime 
and only those with extremely short tidal circulariza- 
tion timescales escape the perturbations of nearby plan- 
ets, we would indeed see a lack of super-eccentric proto- 
hot Jupiters in a sample that lacks young stars. In this 
framework, the observed moderately-eccentric Jupiters 
would not have been produced by HE M bu t by some 
other mechanism, as explored in Section 5.1 Their sur- 



vival indicates that the circularization timescale is not 
extremely short in 3 < P nn ai < 10 days, and therefore 
the argument that planets can only travel along paths 
with very short circularization timescales would only ap- 
ply for hot Jupiters with P < 3 days. 

Finally, a related issue is whether the ratio of 
moderately-eccentric Jupiters to Jupiters with P = P nna i 
is the same in the Kepler and calibration samples (i.e. 
N m od/N P=Pfina ,i ~ iV modi0 /A r p=Pnnai,o), as assumed 
in Equation We would expect the pile-up of hot 

Jupiters to be greater in the older sample (i.e. the Ke- 
pler sample) and therefore that Equation ([6]) might over 
predict the number of super-eccentric Jupiters. However, 
we note that the overall occurrence rate of hot Jupiters 
in the Keple r sample is actually smaller than the in the 
RV-sample ( |Howard et al.|)2012| [Wright et al.p012| , so 
this is unlikely to be a problem in practice. 

5.5. Upper limit on stellar Kozai contribution 

If the observed moderately-eccentric proto-hot Jupiters 
began beyond the ice line with e — > 1 and underwent tidal 
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Fig. 9. — Monte Carlo simulation of tidal evolution of proto-hot Jupiters assuming that all the proto-hot Jupiters were created at once 
(left panel) or that there is a steady current (right panelu). In each simulation, proto-hot Jupiters are drawn from a distribution uniform 
in eccentricity and semi-major axis, extending to 10 AU. Planets that began interior to 1 AU are marked as open, gray circles. We overplot 
tracks of constant angular momentum (dotted lines) corresponding to Pfl na i = 3, 5, 10 days, as well as lines defined by a constant "orbital 
change timescale" (dashed orange lines). In the right panel, the outer orange-dashed line represents a timescale 40 times longer than the 
inner orange-dashed line. The red squares (blue x) were created two (thirteen) times more recently than the black diamonds. 

S12 predicted that if high-eccentricity migration 
(HEM) is the primary channel for producing hot 
Jupiters, the Kepler candidate collection should harbor 
a population of super-eccentric Jupiter-sized planets that 
are in the midst of tidal circularization. We developed 
and performed a procedure to use the publicly-available 
Kepler transit light curves to confirm or rule out this 
prediction and found a paucity of proto-hot Jupiters on 
super-eccentric orbits. Incorporating uncertainties due 
to counting statistics, uncertainties in the light curve fit 
parameters and in the assumed stellar parameters, and 
incompleteness due to the limited observational times- 
pan, we expected to observe 7^3 (Section [3]) but instead 

found only CTtJ (Section j4|). We showed that observa- 
tional effects, including false positives and incomplete- 
ness around noisy stars, are unlikely to explain this dis- 
crepancy. The lack of super-eccentric proto-hot Jupiters 
may indicate that the assumed constant tidal time lag 
approximation — which sets the ratio of super-eccentric 
proto-hot Jupiters to the observed, partially circularized 
moderately-eccentric Jupiters used to compute the pre- 
diction — is incorrect (Section 2.2). However, violation 



circularization while staying on a track of constant angu- 
lar momentum (black arrows, Figure [7]), we would expect 
to see 7l 3 super-eccentric proto-hot Jupiters (Section 3) ; 
the lack of such planets indicates that one or more of the 
alternative pathways in Figure[7] (white, purple, and gray 
arrows) may dominate. These alternative pathways all 
originate from a planetary perturber or disk, rather than 
a stellar perturber. Here we place an upper limit on 
the fraction of hot Jupiters that followed the black ar- 
row channel of HEM, beginning with a super eccentricity 
and moving along a track of constant angular momen- 
tum. Since this is the only pathway open to hot Jupiters 
produced by stellar binary Kozai oscillations, the upper 
limit is also on the fraction of hot Jupiters created by 
stell ar b inaries. We repeat the MCMC procedure in Sec- 
4.2 but update Equation ^ with an additional pa- 



tion 



rameter, representing the fraction of hot Jupiters that 
undergo HEM from super-eccentricities (or, equivalcntly, 
the maximum fraction produced by stellar Kozai): 



iV s 



f*r{e max )N 

mod ■ 



(16) 



In order to fit /* so that N sup agrees with the ob- 
served number of super-eccentric Jupiters -ZV sup , we add 
two terms (for Range 1 and for Range 2) to the log like- 
lihood function (Equation 10 ) of the form: 



AC = N sup ln(7V sup ) - ln(N sup \) - N s 



(17) 



where N sup is the posterior of observed number of super- 
eccentric Jupiters in that Pgnai range (incorporating un- 
certainties in e of the planets in our sample). 

We obtain a posterior distribution for /* that peaks at 
O-lSio'ii' an u PP er limit on the fractional contribution 
from stellar binaries. Therefore, we expect the majority 
of hot Jupiters were created by a planetary perturber 
(or a disk). We note that this result technically only 
applies to hot Jupiters with P > 3 days, because no 
super-eccentric Jupiters with Pgnai < 3 days were ex- 
pected. This limit is consistent with 30% contribution 
rate found by Naoz et al. (2012). 



6. CONCLUSION 



of this assumption could only account for the discrepancy 
if tidal dissipation were actually much stronger at high 
eccentricities along a given angular momentum track. 

In Section [5] we explored a number of dynamical ex- 
planations for the paucity of super-eccentric proto-hot 
Jupiters, relaxing the assumptions by S12 that proto-hot 
Jupiters begin HEM beyond the ice line and that a steady 
current of hot Jupiters is being produced. We found that 
the lack of super-eccentric planets could be explained by 
one of the following scenarios. First, hot Jupiters with 
P > 3 days may be directly implanted interior to the ice 
line, and only those with P < 3 days have undergone 
tidal circularization. This would be the case either if the 
tidal circularization timescale is typically only less than 
a stellar lifetime for Pfj na i < 3 days or if only proto-hot 
Jupiters with fast circularization timescales can manage 
to complete their circularization without being moved by 
a nearby planetary perturber. Second, hot Jupiters with 
P > 3 days may have undergone tidal circularization 
but bypassed the super-eccentric phase of HEM, begin- 
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ning their tidal circularization interior to 0.5 AU with 
0.6 < e < 0.9 rather than beyond the ice line with 
e — > 1. Third, the moderately-eccentric Jupiters used 
to calibrate the prediction may be undergoing Kozai ec- 
centricity oscillations caused by a nearby planetary per- 
turber and we are observing them in the low-eccentricity 
phase, in which they are not currently undergoing tidal 
dissipation. In that case, we would not expect to ob- 
serve super-eccentric Jupiters currently undergoing tidal 
dissipation. 

All these explanations point either to disk migration 
or to secular chaos, planet-planet scattering, or plane- 
tary Kozai (or other yet-to-be-proposed dynamical mech- 
anism) as the dominant channel for hot-Jupiter produc- 
tion, rather than the stellar Kozai mechanism. In Section 
5, we placed an upper limit of 0.15+q'^ on the contribu- 
tion of stel lar Kozai to hot J upiters, consistent with the 
findings of Naoz et al. (2012). Our limit only applies to 
hot Jupiters with orbital periods greater than 3 days, as 
the prediction for super-eccentric Jupiters only applied 
to those ending their HEM journey at 3 < P < 10 days. 

In this paper, we explored S12's prediction for proto- 
hot Jupiters, but they made a similar prediction for a 
population super-ecce ntric binary stars, which they sub- 
sequently discovered ( |Dong et al. 2012). It would not 
be surprising if short-period stars were produced by the 
Kozai mechanism but short-period planets primarily by 
scattering and chaos, which can potentially deliver the 
planets observed interior to the ice line without the plan- 
ets undergoing a super-eccentric phase. The initial con- 
ditions for stellar systems and planetary systems may 
differ in that planetary sy stems are both theo rized (e.g. 
Barnes & Raymond 200 4| and observed (e.g. [Wright et 
al~1|2009| [LTssauer et al.||2(Jll| |Mayor et al.||201l[ ) to form 
packed with many planets, a condition that may often 
lead to scattering and secular chaos. In contrast, stellar 
multiples are typically hierarchical, an optimal setup for 
the Kozai mechanism. 

The lack of super-eccentric proto-hot Jupiters is a new 
piece of evidence that models for making hot Jupiters 
must match, jo ining the distribution of spin orbit mea- 
surcmcnts (e.g. |Fabrycky fc Winn|2009{ |Morton fc John-| 
son|2011a JNaoz et al.|2012p . We recommend that future 
theoretical studies of dynamical models for forming hot 
Jupiters predict: the distribution of timescales for in- 
stabilities that form proto-hot Jupiters, how often the 
high-eccentricity phase of HEM is bypassed, the initial 
conditions along the HEM angular momentum tracks, 
and the expectations for high-eccentricity "failed" hot 
Jupiters that likely have periapses too distant to undergo 
tidal circularization, such as KOI 1474.01 (D12). For the 
brightest Kepler host stars, we recommend measuring 
the spin-orbit alignment of planets in the period valley, 



whose obliquities have presumably not been affected by 
tides. Such measurements could elucidate whether the 
planets in the period valley have a single origin or if there 
are two populations, which might correspond to the cir- 
cular planets and the eccentric planets. Additionally, we 
recommend investigating whether a gas disk could flatten 
and circularize a period valley planet's orbit if the planet 
were scattered there before the gas disk dissipated. 

We recommend that observers strive to better charac- 
terize the eccentricity distribution of the period valley, 
which we argued may be the launching point for the typ- 
ical hot Jupiter's HEM journey. It would be helpful to 
assess if the occurrence rate of eccentric Jupiters in this 
region is - when extrapolated to the 3 < Pfl na i < 10 days 
region - sufficient to launch all the hot Jupiters interior 
to, rather than at or beyond, the ice line. We also recom- 
mend that observers attempt to nail down the period or 
semi-major axis cut-off between the hot Jupiter pile-up 
and the period valley. Finally, although we found that 
it would be more difficult to identify proto-hot Jupiters 
with 0.6 < e < 0.9 using the "photoeccentric effect," it 
could be feasible with more accurate and precise stellar 
parameters. We recommend spectroscopic follow-up of 
KOI host stars for this purpose. 
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APPENDIX 

A. ASSUMPTIONS THAT CANNOT EXPLAIN A LOWER THAN EXPECTED NUMBER OF SUPER-ECCENTRIC 

PROTO-HOT JUPITERS 

These assumptions cannot explain observing fewer than expected super-eccentric proto-hot Jupiters, either because 
a violation would result in more super-eccentric progenitors (1) or because they are unlikely to be violated (2-5). 

1. The evolution of the planet's radius due to tidal inflation is negligible, and no planets are disrupted by tides. 
However, if the planet's radius were to expand over the course of HEM due to tidal inflation, then the tidal 
dissipation rate would be even lower during the earlier stages of HEM, causing planets to spend even longer at 
high eccentricities. Therefore, this effect could only increase the expected number of super-eccentric proto-hot 
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Jupiters. The prevalence of tidal disruption does not affect the S12 prediction, because the prediction is based 
on the survivors. Depending on the timescale of tidal disruption, we may observe additional doomed proto-hot 
Jupiters that will not survive their HEM. 

Angular momentum is not exchanged between the planet and star. If planets were to typically transfer angular 
momentum to stars, we would expect more super-eccentric hot Jupiters than predicted and vice versa. However, 



Penev et al. (2012) argue that stellar tidal dissipation is likely unimportant, because if it were, most hot Jupiters 
would be subsumed by their stars on short timescales. We note that although a star c an add or rem ove angular 
momentum from the planet's orbit as the star rapidly expands on the giant branch (e.g. |Spiegel|2012[ ), we strictly 
restrict our samples to main-sequence stars so we can ignore this effect. We note that the Sun's spin angular 
momentum ranges from about 10% (Pfi na i = 3 days) to 6% (Pfmai = 10 days) the orbital angular momentum of 
a proto-hot Jupiter. 

3. The planet's orbital angular momentum and spin angular momentum are not exchanged. We neglect this effect 
because we assume that the planet maintains a pseudo-synchronous spin throughout its evolution. If the ratio 
of the planet's orbital angular momentum to its spin angular momentum is large, the planet's spin quickly 
(compared to the circularization timescale) evolves to this pseudo-synchronous state, in which the planet's spin 
rate is similar to the orbital frequency at periapse. We expect the ratio of orbital to spin angular momentum is 
indeed typically large, because the planet's distance from the star is very large compared to the planetary radius. 

4. Moderately-eccentric calibration Jupiters (N mo ,i,o) truly have e > 0.2; they are not low- eccentricity planets that 
appear eccentric due to eccentricity bias. Eccentricity bias occurs when noise masquerades as eccentricity. Because 
the eccentricity cannot be negative, it is biased toward higher values. If one decomposes the RV signal caused 
by an eccentric planet into sinusoidal harmonics of the planet's orbital frequency, one finds that the signal due 
to eccentricity i s primarily embedded in the second harmonic and has an amplitude of eK, where K is the RV 
amplitude (e.g. Anglada-Escude et al.|[20T0 ). Eccentricity bias is primarily a concern when eK is near the noise 
level, i.e. for low-mass and/or long-period planets with small K. In contrast, Jupiter-mass planets on short- 
period orbits have large K. For an RV precision of a few m/s and a typical hot Jupiter K ~ 100 m/s, a signal of 
amplitude e x K = 0.2 x 100 m/s = 20 m/s is well above the noise level. Moreover, an even tighter constraint 
on the planet's eccentricity is possible through a joint fit to the RVs and transit light curve, as performed for 
each member of the calibration sample. Therefore we expect that the calibration sample moderately-eccentric 
Jupiters (which have orbital periods ranging from 3-15 days) truly do have e > 0.2. 

5. Only a small fraction of Kepler hot Jupiters are false positives. The expected number of proto-hot Jupiters is 
proportional to the true rate of Kepler hot Jupiters Np—pa na i (Equation p|. For example, if half the Kepler 



hot Jupi t ers wer e fal se positives, the pre dicted number of proto-hot Jupiters should be cut in half. | Morton fc 
Johnson ] (1201 lbh and | Desert et al-1 J2012[ ) find low false-positive rates for Kepler candidates (< 10%). 
et al.| ( |2012[ ) and |Colon et al.| ( |2012[ ) find higher false-positive rates. However, Santerne et al. (20121 focused 
on a population with high a pri ori false-posi t ive pr obabilities because of their V-shaped ligh t curves. Mo r eover 



the false positive discoveries by Colon et al. (2012) were for planets with P < 3 days, which [Colon et aL (2012) 
suggested can be expected from the period distribution of binaries, and we do not include planets with F < 3 days 
in our sample here. We have removed known false-positives from the computations in Section [3j 

B. AVOIDING PROBLEMS DUE TO INCORRECT STELLAR PARAMETERS 
Characterizing the entire planetary eccentricity distribution from transit light curve parameters can be complicated 



by sys tematic errors and uncertainties in the stellar parameters (e.g. |Moorhead et al.|2011 Kane et al.|2012 [ Plavchan 
2012|. Instead, we simply aim to determine whether or not there are light curves for which p c i rc is 



physically 
244/5*. Such a high density 



et al 

unlikely. For example, a planet transiting at periapse with e = 0.95 would have p c i rc 
would be astrophysically implausible based on our knowledge of stellar evolution. 

We consider three potential problems caused by incorrect stellar parameters. First, we could mistakenly identify a 
planet as being highly eccentric even if it had p C i rc ~ pq if we were to underestimate p* as being very low. Second, we 
could miss an eccentric planet if we thought its host star had p* ~ p Q but the true stellar density were much smaller. 
However, we avoid both these problems by restricting our samples to exclude giants. All of the p* we derive for Kepler 
hosts in Section 4 are of order 1. Moreover, by excluding giants from the well-characterized, calibration sample of 
stars with both transits and RV measurements (Figure [T]), we did not make any predictions for super-eccentric planets 
orbiting giants, and therefore cannot miss any. 

The third potential problem is that we used the number of hot Jupiters in the Kepler sample, iVp=pfi n ai, as an 
input for predicting N sup (Equation |6| . If a large fraction of the Kepler hot Jupiters orbits stars that are secretly 
giants but slipped past our stella r parameter cuts, th en we might overpredict the expected number of super-eccentric 
proto-hot Jupiters. As shown by Mann et al. (12012 1 and Dressing et al. (in prep), some of the stars classified as M 
or K dwarfs might be giants. However, Mann et al. ( 2012[ ) find that imposing a cut of \ogg > 4, as we do, helps avoid 
this misclassification. With the cut imposed, 97% of cool stars dimmer than Kp = 14 are dwarves (|M ann et al.|2012j ). 
Among our sample of Kepler hot Jupiters, all the stars with T c g < 5714 K have Kp > 14, so it is very unlikely our 
sample harbors many giants masquerading as M or K dwarves. 
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C. EXCLUDED AND EXCEPTIONAL CANDIDATES 



Based on our fits to the transit light curves, we find that eight candidates are likely eclipsing binaries: KOI-617.01, 
KOI-772.01, KOI-1020.01, KOI-1193.01 KOI-1226.01, KO I -1242 01, KOI-1587.01, and KOI-1773.01. Also, KOI-998.01 



was identified as an eclipsing binary by Dodson- Robinson (2012). We exclude these nine binaries from the rest of this 
work. Finally, we exclude KOI-44.01 because we found, as part of our spectroscopic survey of the brightest KOIs 
(Johnson et al., in prep), that its host star has logg < 4, making it excluded under our stellar cuts. This leaves us 
with 31 planet candidates. The orbital periods of KOI-1477.01 and KOI-1486.01 were incorrectly reported, perhaps 
because of missing transits; their true orbital periods are 169.54 days and 127.2824 days respectively. KOI-211.01 is 
on the list of Kepler eclipsing binaries but this may be an error so we keep it in our sample. 
The transits of KOI-1095.01 were particularly noisy. We applied co-detrending vectors before pro cessing with 



AutoKep and imposed a normal prior on the limb darkening parameters using the routine provided by Sing (2010) 
Typically we can distinguish between a large eccentricity and a large impact parameter for Jupiter-sized planets 
(DJ12), but in this case, we could not tell whether the transit is grazing. To more easily explore parameter space, we 
fit the parameter ln/? c ; rc instead of p c irc- In order to maintain a uniform prior on e and w, we correspondingly use 
the parameter lnp C i rc when performing the MCMC fit to obtain the eccentricity posterior (Section 4.2). Including or 
excluding this candidate only affects the significance of the results by 0.1%. 

KOI-433.02 (Figure 5) has an anomalously low p c i rc , but it has a companion in the system. If KOI-433.02 had 
Pfinai < 10 days, its orbit would cross that of KOI-433.01. Therefore, likely either: a) the candidate has a large 
eccentricity but Pg n ai > 10 days, and is not part of the prediction, or b) the host star's density is overestimated. 
Fortunately, as discussed in Appendix |Bj we have designed our approach so that it is not a problem for our assessment 
of the prediction if some host stars have much lower densities than we estimate. 

As part of our spectroscopic survey, we measured spectroscopic parameters using SpecMatch for KOI-193, KOI-211, 
KOI-686, and KOI-1474, which we use instead of the KIC parameters. For these candidates only, we i mpose a prior 



on th e stellar mass, metallicity, and age from a TRILEGAL (TRIdimcnsional modeL of thE GALaxy; Girardi et al. 



20051 synthetic Kepl e r field population, as described in Section 3.3 of D12. 



The Takeda et al. (20071 evolution models only include stars with > 0.7M Q . For a subset of low-mass host 
stars — KOI 398, KOF333, KOI-806, K OI-855, KOI-918, K OI-1095, KOI-1466, KOI-1477, and KOI-1552 — we use 
the Dartmouth stellar evolution models ( Potter et al.|[2008 l instead, sampling to maintain a uniform prior on stellar 
age, mass, and metallicity. 

D. EFFECT OF MISSING Kepler QUARTERS ON COMPLETENESS 

Some stars in the Kepler sample are missing quarters of data. To assess whether this could account for the lack of 
detected super-eccentric proto-hot Jupiters, we perform the following procedure to estimate C comp (e) in Equation(|7]). 
We tabulate a list of observed quarters for each star in our sample of 50 Kepler hot Jupiter hosts (-/Vp = pg na i). We find 
that 15 stars are missing at least one quarter in Q1-Q8. Because hot Jupiters have short periods and can be detected 
even if quarters are missing, we assume that the missing quarters for the hot Jupiter hosts are roughly representative 
of the entire sample. This approximation allows us to avoid tabulating the missing quarters and performing the 
calculation below for hundreds of thousands of stars. 

Next we calculate C com p(e) for each combination of the following: Pfi na i = 3,5,10 days; P max = 1-5, 2 years; 
0.9 < e < e max in increments of Ae = 10~ 4 ; and each of the 50 lists of observed quarters. For each combination, we 
consider 100 phases for the orbital period P = Pfi n al/(1 — e 2 ) 3 / 2 and calculate C com p(e) as the percentage of phases 
for which two or more transits would be observed, given the missing quarters and P max . C , comp (e) reduces to Equation 
([8| when no quarters are missing. Then we use numerical integration to compute r(e max ) (Equation [7]) and average 
over the 50 lists of observed quarters. This leaves us with an r(e max ) for each Pfi na i and P max . 

We find that missing data cannot account for the lack of super-eccentric proto-hot Jupiters. We repeat the procedure 



of Section 4.2 with the updated r(e max ) and find that our confidence only changes from 98.7% to 98.5% for P max = 2 
years and from 97.6% to 97.2% for P max = 1.5 years. 
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